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Abstract 

Possible origins of the molecular absorption discovered in some ULIRGs are 
investigated, based on a 3-D hydrodynamic model of star-forming interstellar 



' (-h'' gas in a galactic central region. The blue-shifted, warm (~ 200 — 300 K), dense 

(> 10^ cm ^) molecular gas suggested by CO absorption in IRAS 08752+3915 



could be caused by the innermost region of the inhomogeneous inter-stellar 

medium (ISM) around a supermassive black hole. The infrequent observations 

^ . of the dense gas with absorption in ULIRGs and Seyfert 2 galaxies could simply 
O . 

Q^ , suggest that the high-density regions occupy only a very small volume fraction 
00 

CN I of the obscuring material. This is naturally expected if the inhomogeneous 

\6 ' 

O ' structure of the ISM is caused by non-linear development of instabilities. The 

O ■ model predicts a turbulent velocity field in the obscuring material, therefore 
> 

X 

■ for the large enough statistical samples. 



blue- and red-shifted gases should be observable with nearly the same probability 
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dynamics 

1. Introduction 

Active Galactic Nuclei (AGNs) are believed to be highly obscured by dusty dense 
gas, and this contributes to huge radiative energy in some ultra-luminous infrared galax- 
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ies (ULIRG) (Imanishi & Dudley 2000). The central engines of type-2 Seyfert galaxies 
are also likely obscured by optically thick molecular gas, for which a torus-like geometry 
is often assumed. However, detailed structures of the ISM in deep centers of ULIRGs 
and Seyfert galaxies are still unclear. There are some indirect ways to infer their ge- 
ometry and distribution using X-ray spectroscopy (e.g., Levenson, Weaver, & Heckman 
2001), optical/infrared spectral energy distribution (SED) in comparison with model 
SED assuming the geometry of the torus (e.g.. Pier & Krolik 1993; Elitzur 2006; Fritz, 
Franceschini, & Hatziminaoglou 2006), and optical/infrared absorption lines (e.g.. Spoon 
et al. 2004; Imanishi et al. 2006; Imanishi et al. 2007; Levenson et al. 2007). Lutz et 
al. (2004) searched the 4.7 /im fundamental ro-vibrational band of CO in 31 bright lo- 
cal AGNs, but found no clear signature of absorption features, even in a Compton-thick 
Seyfert 2 nucleus, like NGC 1068. Recently, the Spitzer Infrared Spectrograph (IRS) re- 
vealed vibration-rotation absorption bands of gaseous C2H2, HON, and CO2 as well as 
silicate absorption toward deeply obscured (U)LIRG nuclei (Spoon et al. 2004; Lahuis 
et al. 2007; Levenson et al. 2007). The absorption lines of C2H2 and HON suggest the 
presence of warm (Tg ~ 200 — 700 K) and dense (riH > 3 x 10^ cm~^) gas (Lahuis et al. 
2007). They suggest that this gas occupies only a small fraction of the nuclear region (~ 
0.01 pc) near the intrinsic mid-infrared source. The CO absorption features observed in 
ULIRGs are not resolved into individual lines, except IRAS 08572+3915, therefore the 
kinematics and structure of the absorbed material are still an open question. 

Geballe et al. (2006) have found a broad CO absorption line toward IRAS 
08572+3915 NW^ using UKIRT/CGS4. The observations revealed the following fea- 
tures: 1) There are several blue-shifted components with —160 + 25 km s~^ for the high 
J lines (P(6),P(8),...), and -150 + 25 km s~^ and -50 + 25 km s"^ and for the low J 
lines (P(l),i?(l) and R{2)). 2) The hydrogen column density is approximately 1.5 x 10^^ 
cm~^. 3) Mean temperature of the warm blue-shifted component is about 200 K. 

^ Optical classification of this galaxy is LINER (Vcillcux ct al. 1999). Imanishi & Dudley 2000, based 
on infrared spectroscopy, claimed that the AGNs of this galaxy and many other non-Seyfert ULIRGs 
are deeply buried in dusty gas. 
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Shirahata et al. (2007) have also detected the CO absorption in IRAS 08572+3915 
NW using Subaru/IRCS. Their results are consistent with those of Geballe et al. (2006), 
showing a blue-shifted component (-160 km s~^) with the temperature 273 K and cold 
gas at 27 K at the systematic velocity, assuming LTE. They also claim that there is 
a red-shifted component (-1-100 km s~^) seen in high J (J > 4) lines, implying higher 
temperature gas (~ 700 K). Column density of the warm component is 3 x 10^^ cm~^. 

These observations suggest that the absorption features do not simply arise from 
a smooth rotating molecular torus, which has been often postulated to explain the type-1 
and type-2 AGNs. However, it is hard to determine the geometry and internal struc- 
tures of the ISM in the central region of the galaxy from this observational information 
alone. Theoretical models used for SED fitting are phenomenological without kinematical 
information, therefore we cannot compare them with the absorption line observations. 

In this paper, we investigate how the observational features suggested by the 
CO absorption in IRAS-08572+3915 NW can be understood in the context of a three- 
dimensional hydrodynamic model of the ISM around a supermassive black hole (SMBH) 
with nuclear starbursts (Wada & Norman 2002) (hereafter WN02). This is currently a 
unique model of the obscuring material around a SMBH on a several tens pc scale, char- 
acterized by a highly inhomogeneous, multi-phase, and turbulent ISM with a globally 
stable, geometrically thick structure. Three-dimensional radiative transfer calculations 
for this model revealed that CO luminosity distribution is also highly non-uniform (Wada 
& Tomisaka 2005). Although this is not a confirmed theoretical model applicable to the 
obscuring matter in all types of AGNs, it is worth verifying whether the features observed 
in the CO absorption lines in the ULIRG can be explained by the full 3-D hydrodynamic 
models. 
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2. Model and Analysis 

2.1. A Hydrodynamic Model of the ISM around a SMBH 

In the WN02 model, mass, momentum, and energy conservation equations with 
the Poisson equation are numerically solved with energy feedback from supernovae in a 
fixed gravitational potential. A rotating gas disk in a time-independent spherical potential 
is solved by 3-D hydrodynamic code. The mass of the BH is Mbh = lO^M©. We also 
assume a cooling function A(Tg) {JiK <Tg < lO^K) with solar metallicity, heating due to 
photoelectric heating, and energy feedback from SNe. We assume a uniform UV radiation 
field ten times larger than the local UV field. 

The hydrodynamic part of the basic equations is solved by AUSM (Advection 
Upstream Splitting Method) (Liou & Steffen 1993). (See details in Wada (2001); Wada & 
Norman (2001)). We use 256^ x 128 Cartesian grid points covering a 64^ x 32 pc^ around 
the galactic center (i.e. the spatial resolution is 0.25 pc). The Poisson equation is solved 
using the fast Fourier transform and the convolution method. The initial condition is 
an axisymmetric and rotationally supported thin disk (the scale height is 2.5 pc) with a 
uniform radial density profile and a total gas mass of Mg = 5 x 10 ''Mq. Random density 
and temperature fluctuations, which are less than 1 % of the unperturbed values, are 
added to the initial disk. 

Supernova (SN) explosions are assumed to occur at random positions on the disk 
plane. The average SN rate is ~ 0.8 yr'^ The energy of 10 ergs is instantaneously 
injected into a single cell as thermal energy. Thus the three-dimensional evolution of 
blast waves driven by the SNe in an inhomogeneous and non- stationary medium with a 
global rotation is followed explicitly, taking into account the radiative cooling. 

WN02 showed that a globally stable concave "torus", in which the gas is highly 
turbulent and inhomogeneous, is formed (See also figure 6). Figure 1 is a phase-diagram 
of the ISM in the torus at a quasi-steady state. Three dominant phases present: Hot 
gases around Tg ~ 10"^ K caused by SNe, warm gases around Tg ~ 10^ K, and cold, dense 
gases at Tg < 100 - 1000 K. In terms of the observed CO absorption in IRAS 08572-f3915, 
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Fig. 1. Density- Temperature diagram of the ISM around a supcrmassive black liole witli energy 
feedback from supernovae. Tlie open circles represent 500 sampled point in the torus. 

we are interested in kinematics and spatial location of the cold, dense media. 

The size of the present torus model (r ~ 30 pc) is much larger than the size (~ 2 
pc) suggested by near- and mid-infrared high-resolution observations of type-2 Seyfert, 
e.g. NGC 1068 (Jaffe et al. 2004; Wittkowski et al. 2004). This difference actually 
reflects a fact that near /mid-infrared flux is originated from inner a few pc region of the 
dusty material around the AGN. In fact, a recent 3-D radiative transfer model assuming 
clumpy dust tori for NGC 1068 (Honig et al. 2006) suggests that the outer radius of 
the torus that fits SED of NGC 1068 is 56 pc, while their H-, K-, and N-band images 
show their radii are 2.0-2.7 pc. In terms of ULIRGs, Levenson et al. (2007) suggested 
observed deep absorption features in mid-infrared indicate obscuration on scales of a few 
100 pc. Probably the obscuring material in the central region of AGNs and ULIRGs is 
extended from a sub-pc to a sub-kpc scales with various physical/chemical conditions and 
structures. 
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2.2. Analysis of the Hydrodynamic Results 

In order to find the blue- and red-shifted warm, dense gases suggested by the 
CO absorption line in the hydrodynamic model, we sample points with the following 
conditions from the 2000 randomly selected points in a half (i.e. ^ > 0) the "torus". 

1. Number density is n > 10^ cm~^ 

2. Column density from the nuclear source to the observer > 10^^ cm~^ 

3. Temperature is 10 < Tg < 1000 K 
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Fig. 2. Density and temperature of gases wiiicii are candidates for the absorbed gas. 



Figure 2 shows number density and temperature of the selected points. It is clear 
that two components exist: cold gas with ^ 100 K and warm gas with 300 ^ 500 
K. Figure 3 shows the line-of-sight velocity of the selected points, Vios, as a function of 
their temperature. Most of the cold and warm components are distributed in the range of 
l^iosl ^ 200 km s~^. Spatial distribution of the selected points are shown in figures 4 and 
5. The warm components are distributed around 2 ^ r ^ 4 pc from the galactic center, 
and they are located just above the equatorial plane, i.e. 10 ^ 16*1 ^ 30°. Temperature of 
the warm component is higher for larger |^|, i.e. closer to the surface of the "torus". The 
cold dense gases are located 2 ^ r ^ 10 pc and near the disk plane {\9\ ^ 10°). 
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Fig. 3. Same as figure 2, but temperature vs. line-of-sight velocity. 
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Fig. 4. Radial distribution of the selected points. 
In summary, the absorbed material with > 10^^ cm~^, n > 10^ cm~^, |l^os| = 
100 - 200 km s-\ and Tg ^ 300 - 500 K sits at the innermost region of the inhomogeneous 
torus as shown in figure 6 (the hatched regions). For the gas near the equatorial plane, 
the column density is too high to be observed as an absorption line for the continuum IR 
source. Therefore, if the ISM in the central region of ULIRGs, like IRAS 08572+3915, 
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Fig. 5. Same as figure 4, but for temperature vs. position angle 6. = is on the equatorial plane 
of the torus. 

is similar to the structure shown in figure 6, the viewing angle to observe the warm and 
dense g absorption is restricted. If this is the case, we can expect to observe high 

velocity, warm components as well as low velocity, cold components as absorption against 
the background continuum source, whose size is smaller than 1 pc. The high density 
gases in those regions have various velocity components, as shown in figure 3. Thus it is 
not necessary to observe only blue-shifted gas. In fact, Shirahata et al. (2007) also found 
a red-shifted gas component (+100 km s^^) in high- J transition, suggesting gas warmer 
than the blue-shifted component. 

One should note that geometry of the "torus" in the present model depends on 
the supernovae rate, total gas mass, and the mass of a SMBH, since it is determined from 
a balance between energy input and dissipation rates under the effect of gravitational 
potential (WN02). Therefore, the absolute locations (i.e. r and 6) of the warm, dense gas 
around an AGN should depend on these parameters. 
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Fig. 6. (upper panel) Possible location (the hatched region) of the warm gas clumps suggested by 
CO absorption observations. See the lower panel, in which log-scaled gas density is shown at x-y 
and x-z planes (unit of density and scale are Mopc"'^ and 32 pc). 

3. Discussion 

3.1. Volume- Filling Factor of the Warm, Dense Gas 

In the present model of the torus, the high density (> 10® cm~^) gases have both 
red- and blue-shifted velocity for a line-of-sight. Since a volume fraction of these dense 
gases is expected to be very small (see discussion below), and the background continuum 
is emitted from a small region in the vicinity of the central engine, the chance to observe 
a dense blob with a particular velocity component by an absorption line should be very 
small. For example, as shown in figure 3, the blue-shifted component around —200 km 
s~^ exists in about 1/500 of the total volume of the torus whose radius is about 30 
pc. Therefore, one possible explanation about the observed blue-shifted component is 
that the line-of-sight toward the AGN happens to graze the high-density clump which 
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is approaching us. If this is the case, since the motion in the torus is highly turbulent 
(WN02), we have nearly equal chance of observing red-shifted warm gas, too. The red- 
shifted warm gas found by Shirahata et al. (2007) may be explained by this picture, but 
observed samples are still too small to further compare with the model. 

Recently Wada & Norman (2007) proposed a simple statistical theory of the struc- 
ture of the inhomogeneous gas disk in galactic disks, which is caused by non-linear devel- 
opment of gravitational and thermal instabilities. Using three-dimensional hydrodynamic 
simulations, they show that the probability distribution function (PDF) of density P{p) in 
a globally stable, inhomogeneous ISM is well represented by a single log-normal function 
over a wide density range, that is 



P{p)dp = ^=^exp 



2a> ^InP, 

where po is the characteristic density and a is the dispersion. The simulations show that 
the dispersion of the log-normal PDF is larger for more massive systems. Using the PDF, 
it is straightforward to calculate a volume fraction of high density gas (/c) above a given 
critical density (pc): 

ln(Pc/po) -f^ 



«p„.) = i(l-Erf 
The dispersion a is then related to the average density, p, i.e. 



(2) 



It is also suggested that characteristic density po is not very sensitive for changing the 
total gas mass. Suppose the column density toward the continuum source is 10^^ cm~^ 
and the size of the 'torus' is 30 pc, the average density would be p ~ 10^ cm~'^. Therefore, 
using equation (3), we have a ~ 3.0 for po = 1 cm~'^. This leads to /c = 9.5 x 10"^ and 
fc = 5.3 X 10"^ for Pc = 10^ cm~^ and 10^ cm~^, respectively. With this small volume-filling 
factor of the high-density gas in an obscuring material, it might be reasonable to expect 
there is a small chance of detecting a high-density clump with a particular velocity by 
absorption lines. 
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3.2. Difference between ULIRG and Seyfertl,2 

Hao et al. (2007) observed 196 AGNs and ULIRGs using 5j»zfeer/IRS, and found 
that quasars are characterized by sihcate features in emission and Seyfert2s are dominated 
by weak sihcate absorption, while ULIRGs are characterized by strong sihcate absorption. 
Levenson et al. (2007) suggest that this difference in near-infrared spectra reflects a differ- 
ent geometry of the dusty material. They claim that clouds dominate the obscuration of 
Seyfert Is and 2s, while on the other hand, central sources of ULIRGs are obscured mainly 
by smooth, geometrically and optically thick material with a steep temperature gradient. 
According to Lahuis et al. (2007), A^H2 fo^^ the nucleus of IRAS 08572+3915 estimated by 
infrared spectra obtained by Spitzer is N^^ ^ 1.5 x 10^^ cm~^. If the scale of the obscuring 
material is about 30 (100) pc, the average number density is n ^ 1.5 (0.5) x 10^ cm~^. 
Using the log-normal density PDF in §3.1, we can estimate that the volume filling factor 
of those gases is about 20 (30) %, suggesting that the dusty gas that contributes to the 
silicate absorption is more smoothly distributed in the obscuring material than the warm, 
dense gas found by CO absorption. 

Our result presented here is that at least the CO absorbed feature of IRAS 
08572+3915 is not inconsistent with a model characterizing a dense clumpy medium 
around an AGN. Since the velocity field of the ISM is highly turbulent in the model, if 
individual high density clumps are not resolved, a broad absorption feature originating in 
dense gases with various velocity components can be expected in the obscuring material. 
If this is the case, the reason why CO absorption is not observed in Seyfert-2 (Lutz et al. 
2004) can be understood by the low filling factor of the warm, dense gas in the obscuring 
material. As discussed in §3.1, if average gas density of the obscuring material around 
AGNs is smaller in Seyfert galaxies than in ULIRGs, a fraction of high-density clumps be- 
comes small. This seems to be inconsistent with the result by Levenson et al. (2007), but 
probably dusty "clouds" in the obscuring material responsible for the near/mid-infrared 
SED are much more diffuse than those causing the CO absorption (see discussion above). 
The ISM around AGNs in ULIRGs and Seyfert galaxies are more or less inhomogeneous. 
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and if so its density PDF is expected to be approximately log-normal. Therefore, it is not 
surprising, if the "smooth" distribution of dusty gas that is preferable to the SED fitting 
would include higher-density clumps with a smaller volume-filling factor. 

3.3. Limitations of the Current Model and Other Possibilities 

In the present hydrodynamic model, heating due to the radiation from the AGN 
to dust is not taken into account. Temperatures of optically thick clouds at a distance of 
a few pc from an AGN with L = IO^^^q would be 300-1200 K (Elitzur 2006). Therefore 
the warm clouds near the surface of the clumpy "torus" might be warmer than 300 — 500 
K, which is a result of the present model. Effects of X-rays are not included in the model 
either, which would be especially important for chemistry of molecules (Maloney et al. 
1996; Meijerink et al. 2007). For highly inhomogeneous media like those seen in figure 6, 
hard X-rays could affect not only chemistry of the gas in the vicinity of the central source, 
but also the clouds inside the "torus" , because the column density and the ionization rate 
are not simply related to the distance from the central source (WN02). The effects of 
X-rays on the abundance and line intensities of CO and other molecules are investigated 
using the present hydrodynamic model and 3-D non-LTE radiative transfer calculations 
(Yamada, Wada, & Tomisaka, in preparation). 

As an alternative interpretation of the blue-shifted warm gas could be originated in 
outflows from an accretion disk around a supermassive black hole. Using two-dimensional 
radiation-hydrodynamic simulations, Ohsuga (2006) found that an outflow is driven by 
radiation force due to a luminous accretion disk. However, the outflow velocity is ex- 
tremely large (~ 0.1c), and its temperature is ~ 10^ K, therefore the outflow from the 
accretion disk itself does not explain the warm, blue-shifted gases discussed here. It might 
be possible however that, the hot outflow gas propagates outward, and it then cools and 
forms warm, dense cloudlets by the thermal instability at a few pc from the accretion 
disk. Unfortunately long-term evolution of the outflow from the accretion disk and its 
interaction with the ISM on several tens pc scale is still an open question. 
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4. Conclusion 



Based on the analysis presented here, we could conclude that the blue-shifted, 
warm, dense molecular gas in IRAS 08752+3915 NW suggested by recent high-resolution 
observations by Subaru/IRCS can be caused by high-density clumps in the inhomogeous 
inter stellar matter around a supermassive black hole. The small possibility of observing 
the CO absorption in ULIRGs and Seyfert 2 galaxies could be simply interpreted as indi- 
cating that the high density regions occupy a very small volume fraction (< 10^'^ — 10"^) 
in the obscuring material with a turbulent velocity field. This is naturally expected if the 
inhomogeneous material is caused by non-linear development of instabilities and random 
processes (e.g. stochastic explosions of supernovae and interaction with inhomogeneous 
ISM) in a globally quasi-stable system. 

However, this model is still a presumption until the absorbed features in near in- 
frared spectra of many ULIRGs and Seyfert 2 galaxies are resolved and details of the 
kinematics of the obscuring material are clarified by future observations. It is also neces- 
sary to improve numerical models of the ISM in the central region of active galaxies by 
taking into account more realistic treatment of chemical reactions of molecular gas and 
radiative transfer in terms of AGNs and star forming regions. 

The author is grateful to Takao Nakagawa, Mai Shirahata, and Masa Imanishi for 
fruitful discussions. The anonymous referee's comments were also helpful. Numerical 
computations were carried out on a Fujitsu VPP5000 at NAOJ. KW is supported by 
Grant-in-Aids for Scientific Research [no. 16204012] of JSPS. 
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